Lattice study of exotic S — +1 baryon 
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We propose 5=+l baryon interpolating operators, which are based on an exotic description of the 
antidecuplet baryon like diquark-diquark-antiquark. By using one of the new operators, the mass 
spectrum of the spin-1/2 pentaquark states is calculated in quenched lattice QCD at (3 — 6/g 2 = 6.2 
on a 32 3 x 48 lattice. It is found that the J p assignment of the lowest O (uudds) state is most likely 
(l/2)~. We also calculate the mass of the charm analog of the O and find that the Q c (uuddc) state 
lies much higher than the DN threshold, in contrast to several model predictions. 
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Recently, LEPS collaboration at Spring-8 has observed 
a very narrow resonance + (154O) in the K~ missing- 
mass spectrum of the jn — > nK + K~ reaction on 12 C 
A remarkable observation is its strangeness is S=+l, 
which means that the observed resonance must contain a 
strange antiquark. Thus, the O + (1540) should be an ex- 
otic baryon state with the minimal quark content uudds. 
This discovery is subsequently confirmed in different re- 
actions by several other collaborations 0- It should be 
noted, however, that the experimental evidence for the 
+ (154O) is not very solid yet since there are a similar 
number of negative results to be reported Q . 

Theoretically the existence of such a state was pre- 
dicted long time ago by the Skyrme model 0]. However, 
the prediction closest in mass and width with the ex- 
periments was made by Diakonov, Petrov and Polyakov 
using a chiral-soliton model 0- They predicted that it 
should be a narrow resonance and stressed that it can be 
detected by experiment because of its narrow width. In 
a general group theoretical argument with flavor SU(3), 
S=+l pentaquark state should be a member of antidecu- 
plet or higher dimensional representation such as 27-plet 
or 35-plet. Both the Skyrme model and the chiral-soliton 
model predict that the lowest 5=+l state appears in 
the antidecuplet, 1=0, and its spin and parity should be 
Experimentally, spin, parity and isospin of 
the + (154O) are not determined yet. After the discovery 
of the + (154O), many model studies for the pentaquark 
state are made with different spin, parity and isospin. 

Lattice QCD in principle can determine these quantum 
numbers of the + (154O), independent of such arbitrary 
model assumptions or the experiments. We stress that 
there is substantial progress in lattice study of excited 
baryons recently 0- Especially, the negative parity nu- 
cleon TV* (1535), which lies close to the 6+ (1540), has be- 
come an established state in quenched lattice QCD |a,u|- 
Here we report that quenched lattice QCD is capable of 
studying the 6 + (1540) as well. 

Indeed, it is not so easy to deal with the qqqqq state 
rather than usual baryons (qqq) and mesons (qq) in lat- 
tice QCD. The qqqqq state can be decomposed into a pair 
of color singlet states as qqq and qq, in other words, it 



can decay into two-hadron states even in the quenched 
approximation. For instance, one can start a study with 
a simple minded local operator for the + (1540), which is 
constructed from the product of a neutron operator and a 
K + operator such as = s a bc{d^,C^Ub)d c {s e , y5U e ) , The 
two-point correlation function composed of this operator, 
in general, couples not only to the state (single hadron) 
but also to the two-hadron states such as an interacting 
KN system 0, Even worse, when the mass of the 
qqqqq state is higher than the threshold of the hadronic 
two-body system, the two-point function should be domi- 
nated by the two-hadron states. Thus, a specific operator 
with as little overlap with the hadronic two-body states 
as possible is desired in order to identify the signal of the 
pentaquark state in lattice QCD. 

For this purpose, we propose some local interpolating 
operators of antidecuplet baryons based on an exotic de- 
scription as diquark-diquark-antiquark. There are basi- 
cally two choices as 3 C £g> 3 C or 3 r ® 6 r to construct a color 
triplet diquark-diquark cluster |10Llll| . We adopt the for- 
mer for a rather simple construction of diquark-diquark- 
antiquark. We first introduce the flavor antitriplet (3/) 
and color antitriplet (3 C ) diquark field 

\ x ) = \£i]k£abcq^ b {x)CTq k ^ c (x) (1) 



where C is the charge conjugation matrix, abc the color 
indices, and ijk the flavor indices. T is any of the six- 
teen Dirac 7-matrices. Accounting for both color and 
flavor antisymmetries, possible Ts are restricted within 
1, 75 and 757 p which satisfy the relation (CT) T = — CT. 
Otherwise, the above defined diquark operator is identi- 
cally zero. Hence, three types of flavor 3/ and color 3 C 
diquark; scalar (75), pseudoscalar (1) and vector (757^) 
diquarks are allowed hj. The color singlet state can be 
constructed by the color antisymmetric parts of diquark- 
diquark (3 C (g> 3 c )antisym = 3 C with an antiquark (3 C ). In 
terms of flavor, 3/<g)3/<g>3/ = 1/ ® 8/ ffi 8/ ffilO/. Mani- 
festly, in this descri ptio n, the 5=+l state belongs to the 
flavor antidecuplet [13J. Automatically, the 5=+l state 
should have isospin zero. Then, the interpolating opera- 
tor of the (uudds) is obtained as 
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FIG. 1: Examples of the effective mass for (a) the positive parity Q(uudds) state and (b) the negative parity 0(uudds) state 
at k = 0.1506 for up and down quarks and the strange quark k s — 0.1515. The solid lines represent fitted mass and their 
statistical errors. The dashed lines correspond to the total energy of the non-interacting KN state with the smallest possible 
nonzero lattice momentum for (a) and with zero momentum for (b). 



for r / T'. The form Cs T for the strange antiquark 
field is responsible for giving the proper transformation 
properties of the resulting pentaquark operator under 
parity and Lorentz transformations Note that be- 

cause of the color antisymmetry, the combination of the 
same types of diquark is not allowed. Consequently, we 
have three different types of exotic $=+1 baryon oper- 
ators through the combination of two different types of 
diquarks, which have different spin-parity |l2| : 

G>\ = eabc£aef£b g h(u^Cd f )(u^C-f 5 d h )Cs^, (3) 
= £ abc e ae f£bgh(UeCj 5 df)(UgC7 5 jpd h )CSc,(4:) 
Bl )A( = ea& c ea e /e6 S / i ("fC'd / )(uJC757 p d/ l )CsJ (5) 

where the subscript "+ (— )" refers to positive (negative) 
parity since these operators transform as VQ±(x, t)V' = 
±-f4<d±(—x,t) (for /1 = 1,2,3) under parity. The first 
operator of Eq. (J3J) is proposed for QCD sum rules in a 
recent paper |14| independently. 

In this description, the operator of exotic S3/2 (ssddu 
or uussd) states, which are members of the antidecuplet, 
can be treated by interchanging u and s or d and s in 
the above operators. If a strange antiquark is simply 
replaced by a charm antiquark, the proposed pentaquark 
operators can be regarded as the anti-charmed analog of 
the isosinglet pentaquark state, Q c (uuddc). 

Recall that any of local type baryon operators can cou- 
ple to both positive- and negative-parity states since the 
parity assignment of an operator is switched by multiply- 
ing the left hand side of the operator by 75 . The desired 
parity state is obtained by choosing the appropriate pro- 
jection operator, 1 ± 74, on the two-point function G(t) 
and direction of propagation in time. Details of the par- 
ity projection are described in Ref. 7j. We emphasize 
that the second and third operators, Eqs. 10} and 10, 
can couple to both spin-1/2 and spin-3/2 states. By using 
them, it is possible to study the spin-orbit partner of the 
spin-1/2 O state, whose presence contradicts the Skyrme 



picture of the 6 [ll|. However, we will not pursue this 
direction in this article. We utilize only the first operator 
of Eq. ((2J), which couples only to a spin-1/2 state. Un- 
der the assumption of the highly correlated diquarks, we 
simply omit a quark-exchange diagram between diquark 
pairs contributing to the full two-point function in the 
following numerical simulations. 

We generate quenched QCD configurations on a lattice 
L 3 xT = 32 3 x 48 with the standard single-plaquette Wil- 
son action at (3 = 6/g 2 = 6.2 (a" 1 = 2.9 GeV). The spa- 
tial lattice size corresponds to La ~ 2.2fm, which may be 
marginal for treating the ground state of baryons with- 
out large finite volume effect. Our results are analyzed 
on 240 configurations. The light-quark propagators are 
computed using the Wilson fermions at four values of the 
hopping parameter k = {0.1520,0.1506,0.1497,0.1489}, 
which cover the range M n /M p = 0.68-0.90. n s = 0.1515 
and k c = 0.1360 are reserved for the strange and charm 
masses, which arc determined by approximately repro- 
ducing masses of 0(1020) and J/*(3097). We calculate 
a simple point-point quark propagator with a source lo- 
cation at t src = 6. To perform precise parity projection, 
we construct forward propagating quarks by taking the 
appropriate linear combination of propagators with pe- 
riodic and anti-periodic boundary conditions in the time 
direction. This procedure yields a forward in time prop- 
agation in the time slice range < t < T — t src . 

In this calculation, the strange (charm) quark mass is 
fixed at n s (k c ) and the up and down quark masses are 
varied from M„ w 1.0 GeV (k = 0.1489) to M„ « 0.6 
GeV (k — 0.1520). Then, we perform the extrapolation 
to the chiral limit using five different n values. 

We first calculate the effective masses M c g(t) = 
ln{G(t) / 'G(t + 1)} for both parity states of the spin-1/2 
<d(uudds). For example, Figs.Hshow effective masses for 
the positive parity channel and the negative parity chan- 
nel at k — 0.1506 for up and down quarks with the fixed 
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FIG. 2: An example of the effective mass for the negative 
parity state of the charm analog ® c (uuddc) at K = 0.1506 for 
up and down quarks and the charm quark k c — 0.1360. 



strange quark. Statistical uncertainties in both figures 
are estimated by a single elimination jack-knife method. 

In Fig. H (a), the effective mass plot for the positive 
parity state shows a very short plateau albeit with large 
statistical errors. This plateau terminates at t w 17 and 
then the rather noisy signals appear after t = 18 and 
become reduced around the KN threshold. We remark 
that the positive parity state can decay into the KN 
state in a P-wave where the two hadrons should have a 
nonzero momentum. However, all momenta are quan- 
tized as p n = 2irri/L on a system of finite volume. The 
KN threshold is defined as the total energy of the non- 
interacting KN state with the smallest nonzero momen- 
tum |pmin| = 2tt / L in lattice units. Here, we stress the 
following two points. First, there is no clear signal for 
the KN state to be observed in the effective mass plot. It 
means that our proposed interpolating operator couples 
weakly to the KN scattering state. Secondly, our ob- 
served plateau in Fig. (a) is considerably higher than 
the KN threshold. While the observed asymptotic state 
can be identified as a pentaquark (single hadron) state, 
our results seem to give no indication of the O + (1540) 
state in the positive parity channel. 

In the negative parity channel, the gross feature is sim- 
ilar to the case of the positive parity. Fig. (b) shows 
that a clear plateau appears in the range 15 < t < 20. 
The relatively noisy signals appear around the S-wave 
KN threshold after t = 21 and continue toward the max- 
imum time slice t = T — t SIC for the forwarding propa- 
gation. The errors after t = 21 are probably underesti- 
mated. The correlators for the heavier mass state in this 
euclidean time region usually have many orders of mag- 
nitude deviation and the distribution is non-Gaussian. 
Therefore, the signals after t = 21 are inconclusive. 

We perform a covariant single exponential fit [22j to 
the two-point function in the plateau region 15 < t < 20, 
where the respective x 2 is indeed most favorable. The 
estimated mass is clearly higher than the KN thresh- 



old, which is evaluated as the total energy of the non- 
interacting KN state with zero momentum. The exci- 
tation energy of the observed asymptotic state from the 
KN threshold is roughly consistent with the experimen- 
tal value. Although, without a finite volume analysis, it 
cannot be excluded that the observed plateau stems from 
only a mixture of the KN scattering states; we may re- 
gard it as a pentaquark state with a mass close to the 
experimental value of the + (154O). 

As a strange antiquark is simply replaced by a charm 
antiquark, we can explore the anti-charmed pentaquark 
Q c (uuddc) as well. A similar identification for the C 
state can be made in the negative parity channel. The ef- 
fective mass plot (Fig. |2| shows that a plateau, which ter- 
minates at t ~ 21, is much higher than the DN threshold. 
The relatively noisy signals appear around this threshold 
after t « 21. The observed asymptotic state is identified 
with a pentaquark (single hadron) state similarly. 

In Fig. |3| we show the mass spectrum of the Q(uudds) 
states with the positive parity (open squares) and the 
negative parity (open circles) as functions of the pion 
mass squared. Mass estimates are obtained from co- 
variant single exponential fits in the appropriate fit- 
ting range. All fits have a confidence level larger than 
0.3 and x 2 /Ndf < 1-2. It is evident that the low- 
est state of the isosinglet $=+1 baryons has the nega- 
tive parity. We evaluate the mass of the Q(uudds) with 
both parities in the chiral limit. A simple linear fit for 
all five values in Fig. yields Men/ 2 -)=0.62 (3) and 
Mq( 1 / 2 +)=1-00 (5) in lattice units. If we use the scale set 
by r from Ref. QjJ, we obtain M e(1/2 - } = 1.84(8) GeV 
and Me(i/2+) = 2.94(13) GeV. It is worth quoting other 
related hadron masses. The chiral extrapolated values for 
the kaon, the nucleon and the N* state are Mk = 0.53(1) 
GeV, M N = 1.06(2) GeV and M N * = 1.76(5) GeV in this 
calculation. 

Our obtained 0(l/2~) mass is slightly overestimated 
in comparison to the experimental value of the O + (1540), 
but comparable to our observed N* mass, which is also 
overestimated. Needless to say, the evaluated values 
should not be taken too seriously since they do not in- 
clude any systematic errors. Such a precise quantitative 
prediction of hadron masses is not the purpose of the 
present paper. Rather, we emphasize that, our results 
strongly indicate the J p assignment of the + (154O) is 
most likely (l/2)~. This conclusion is consistent with 
that of a recent lattice study E3 (if one corrects the par- 
ity assignment of their operator |l7|) and that of QCD 
sum rules approach 

El- 
Results for the lowest-lying spin-1/2 O c state, which 
has the negative parity, are also included in Fig. The 
8 C state lies much higher than the DN threshold in con- 
tradiction with several model predictions 0, 0, l2oT |. 
The chiral extrapolated value of the 6 C mass is 3.45(7) 
GeV, which is about 500 MeV above the DN threshold 
(Mrj=1.89(l) GeV) in our calculation. This indicates 



4 



3.0 



2.5 - 



2.0 



a 

£ 1.5 

US 



1.0 



0.5 



0.0 







@ c (l/2-) 


- 






- 


.33-- 


@(i/2+) 








- ■#■-"" 




....®'"" 

0(1/2-) 


I 


..-•"®' 

i i 


i i 


0.00 


0.05 0.10 

(aMjt) 


0.15 0.20 

2 



FIG. 3: Masses of the spin-1/2 Q(uudds) states with both 
positive parity (open squares) and negative parity (open cir- 
cles) as functions of pion mass squared in lattice units. The 
charm analog Q c (uuddc) state (open diamonds) is also plot- 
ted. Horizontal short bar represents the KN(DN) threshold 
estimated by Afjv + Mk(Mn + Mn) in the chiral limit. 



that the anti-charmed pentaquark O c is not to be ex- 
pected as a bound state. 

We have calculated the mass spectrum of the 5=+l ex- 
otic baryon, <d(uudds), and the charm analog Q c (uuddc) 
in quenched lattice QCD. To circumvent the contam- 
ination from hadronic two-body states, we formulated 
the antidecuplet baryon interpolating operators using an 
exotic description like diquark-diquark-antiquark. Our 
lattice simulations seem to give no indication of a pen- 
taquark in the positive parity channel to be identified 
with the + (154O). In contrast the simulations in the 
negative parity channel can easily accommodate a pen- 
taquark with a mass close to the experimental value. Al- 
though more detailed lattice study would be desirable to 
clarify the significance of this observation, the present 
lattice study favors spin-parity (l/2)~ for the O + (1540). 
We have also found that the lowest spin-1/2 6 C state, 
which has the negative parity, lies much higher than 
the DN threshold, in contrast to several model predic- 
tions dill El. 

To establish the parity of the + (154O), more exten- 
sive lattice study is required. Especially, a finite volume 
analysis is necessary to disentangle the pentaquark sig- 
nal from a mixture of the KN scattering states. It is 
also important to explore the chiral limit. This calcu- 
lation was performed using relatively heavy quark mass 
so that one may worry about a level switching between 
both parity states toward the chiral limit as observed in 
the case of excited baryons 0, 0] . We remark that a 
study for the non-diagonal correlation between our pen- 



taquark operators and a standard two-hadron operator 
should shed light on the structure of the very narrow res- 
onance O + (1540). The possible spin-orbit partner of the 
6 state is also accessible by using two of our proposed 
operators. We plan to further develop the present cal- 
culation to involve more systematic analysis and more 
detailed discussion. 
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A simple exponential might not be an appropriate func- 
tional form for the decaying state. Recall that the pen- 
taquark state can decay into two-hadron states even 
in the quenched approximation. Strictly speaking, we 
should take the decay width into account in the fitting 
form. However, nobody knows an appropriate analytic 
form for the two-point function of an unstable state in 
finite volume on the lattice. 



